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ABSTRACT 



An integrating mierothrust balance was utilized to 
evaluate the impulse and thrust performance characteristics 
of a rocket system which employs a new concept for pro- 
pulsion; that is, it effects a controlled thrust from the sub- 
limation of a solid propellant. The propulsion package 
tested was designed to provide thrust in a control system of 
a spin stabilized satellite. 
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IMPULSE AND THRUST TEST OF A 
SUBLIMATING-SOLID MICROPROPULSION SYSTEM 

by 

Richard W. Forsythe 
Goddard Space Flight Center 


INTRODUCTION 

This report describes the initial tests conducted by Goddard Space Flight Center on a micro- 
thrust attitude control rocket system which employs a new concept; that is, it utilizes a sublimating 
solid propellant as a means for propulsion. The system tested was designed and fabricated on NASA 
contract NAS5-3153 by Rocket Research Corporation of Seattle, Washington. The concept of the 
propulsion developed as a result of in-house R&D work by the same contractor. 

There are a large number of applications for very low thrust systems in satellite missions. 
Until now, the only operational hardware systems available utilized cold gas nitrogen as the pro- 
pellant. Cold gas nitrogen has the major disadvantages of (a) high weight due to low density of the 
propellant and high pressure of the fuel tank and (b) reduced reliability of the system due to 
multiplicity of high pressure moving parts such as pressure regulators and control valves. 

The sublimating solid microrocket is an improved cold-gas reaction jet system, which main- 
tains cold gas simplicity while reducing system component weight and increasing reliability. The 
sub limating ; solid microrocket has a theoretical vacuum specific impulse of 85 seconds at 70° F when 
expanded through a nozzle with a 40 to 1 area ratio. Physically, the system consists of a thin wall 
propellant tank, a sublimating solid propellant in crystal form, a filter, a propellant valve, 
and a nozzle. The resulting control rocket is one of the simplest systems possible. When 
the valve is opened, the vapor escapes through the nozzle to produce thrust. The solid sublimes 
to replace the vapor which has been released. The use of the sublimating solid propellant does 
not involve ignition, combustion, or high temperature and high pressure. Therefore, no reaction 
chambers, pressure regulators, or high pressure storage tanks are required. Since the propellant 
is stored in solid form and the vapor pressure of the propellant is low, the propellant tank weight 
for the sublimating solid microrocket is very low, in direct contrast to the cold-gas nitrogen 
system in which tank weight often exceeds propellant weight. Low system operating pressure also 
reduces leakage problems and results in high overall reliability for the system. 

TEST SYSTEM DESCRIPTION 

The NASA specifications for the propulsion system were selected such that the operating 
characteristics met the requirements of a spinning- satellite-pointing controls mechanism in all 
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respects except weight. The components were designed for laboratory testing only and no attempt 
was made to minimize weight. The sublimating- solid propulsion system tested is shown assembled 
in Figure 1. 

The total weight of the laboratory system was 2.67 pounds with approximately 1.6 pounds of 
fuel. The valve utilized was an axial solenoid poppet type valve. 

The fuel utilized was selected considering the following criteria: 

1. The propellant material must possess a significant vapor pressure (in relation to the desired 
thrust level). 

2. The melting or decomposition points cannot be too low. 

3. The molecular weight of the exhaust gases should be as low as possible to achieve a high 
specific impulse. 

4. The propellant density should be as high as possible to reduce storage area and weight. 

5. The heat of sublimation of the solid fuel should be as low as possible to allow higher 
impulse bits. 

From a number of possible fuels, ammonium hydrogen sulfide (NH 4 HS) was selected as the fuel. 
Ammonium hydrogen sulfide has a heat of sublimation of 782 Btu/lbm and a density of approximately 
0.03 lbm/in 3 . Upon sublimation it decomposes into ammonia (NH 3 ) and hydrogen sulfide (H 2 S) with 
an average molecular weight of 25.6. The gases have a distinct and pungent odor; although high 
concentrations are toxic, due to their low threshold of olefactory sensitivity, they are easily de- 
tected long before the concentrations are sufficient to be dangerous. Figure 2 shows the equilib- 
rium conditions for ammonium hydrogen sulfide on a pressure versus temperature curve. 



Figure 1— Assembled subl imating-sol id 
micropropulsion system. 


Table 1 


System Specifications. 


Parameter 

Specification 

Total Thrust (2 nozzles) 

5 x 10~ 2 lbf (nominal) 

Specific Impulse 

83 lbf-sec/lbm (theoretical) 

Inlet Nozzle Pressure 

1.80 psia (nominal) 

Tank Pressure 

6.35 psia (nominal) 

Total Mass Flow 

6 . 04 x 1 0" 2 lbm/ sec (nominal ) 

Thrust Coefficient 

1.80 (theoretical) 

Throat Area 

7.85 x 10" 3 in 2 

Area Ratio 

40 


2 



SUBLIMATION PRESSURE 



Figure 3 — Fuel tank for sublimating-solid micropropul- 
sion system and its styrofoam tank container. 


Figure 2—Subl imation pressure versus temperature. 

TEST PLAN 




The purpose of the test was to evaluate the 
thrust and impulse delivered from the pro- 
pulsion system which had been stored with its 
fuel in a refrigerator system at 34 F for 8 
months. To eliminate the extraneous thermal 
environment around the test system during the 
tests, the fuel tank and valve were essentially 
styrofoam container. Figures 3 and 4 illustra 


Figure 4 — Styrofoam container completely 
enclosing fuel tank. 

thermally isolated by mounting them inside a 
e the test arrangement. 


To achieve the desired temperatures during the test, heater strip wiring was wrapped directly 
around the fuel tank and valve and supplied with power upon command external to the vacuum 
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chamber. Thermocouples fed into the 
styrofoam box system to allow tempera- 
ture monitoring. The valve temperature 
was kept above the fuel temperature at 
all times to insure against any gas re- 
condensation at the valve. 

To accurately measure the thrust 
and impulse, the propulsion system was 
mounted on an Integrating Microthrust 
Balance, (IMB) in a vacuum chamber 
(Figure 5). The IMB, which utilizes a 
soft suspension, is basically a linear dis- 
placement pendulum designed and instru- 
mented to measure momentum imparted 
to it by a single pulse bit. 

The IMB is calibrated by instantane- 
ously applying a known momentum to the 
pendulum; because the propulsion system is not an instantaneous application of total momentum, 
a slight error in readout develops with the amount directly related to the impulse time. Analysis, 
however, allows the development of correction factors to compensate for the impulse time and also 
drag of the pendulum due to a damping system. All data shown has been adjusted slightly by a 
correction curve. Appendix A gives a complete description of the IMB and the analysis which de- 
velops the correction curve. 

Once the test system and equipment was set up as shown in Figure 5, the vacuum chamber was 
pumped down to 2 x 10' 2 mm of Hg at which pressure three tests were conducted to evaluate the 
thrust and impulse of the system as follows: 

Test 1 — The propulsion system valve was energized and the system allowed to deliver a steady 
state thrust for 18.33 seconds. 

Test 2 — Power to the solenoid valve was de- 
livered by a timing circuit exter- 
nal to the chamber which allowed the 
system to pulse for 250 milliseconds 
every 3 seconds for 40 pulses. 


Test 3 — By the timing circuit, the system was 
pulsed with varying pulse bits 250, 
200, 150, 100 and 77 milliseconds in 
duration. 



Figure 6— Location of thermocouples on fuel tank. 
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Two thermocouples at the positions indicated in Figure 6 monitored the temperature during 
all tests. 

TEST RESULTS AND ANALYSIS 

The following test results and analysis deal with the impulse and thrust data after it was re- 
duced and adjusted by the IMB procedure discussed in Appendix A. 

Test 1 

A power supply external to the chamber .supplied 28.0 volts to the solenoid valve for 18.33 
seconds. This allowed the propulsion system to deliver a steady state thrust which displaced the 
IMB pendulum. The test record of displacement versus time is shown in Figure 7. 



Figure 7— Plot of relative displacement versus time for steady-state thrust. 


Extrapolating the steady- state thrust displacement back to time = 0 seconds, an initial rel- 
ative displacement of 2.28 is obtained. The final relative displacement is 1.40. 

The equation of motion describing the thrust balance motion is 

mx + cx + kx = F(t) . (1) 

Under the influence of a steady state forcing function independent of time where the damping, c, is 
appreciable, the resultant displacement is described by 

kx = F when x and x approach 0 . (2) 

Utilizing this fact and the dynamics of the IMB (see Appendix A for deviation of thrust levels) the 
thrust levels at thrust initiation and termination were obtained. 
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Table 2 presents the Test 1 data. The terminal temperatures shown were the lowest values 
of the temperatures and were reached at approximately T + 125 seconds as indicated. 

Table 2 


Test 1 Data. 


Time (seconds) 

Relative Displacement 

Thrust, 2 nozzles, (lbf) 

Temp No. 1 (°F) 

Temp No. 2 (°F) 

0 

2.28 

0.068 

73.4 

87.0 

18.33 

1.40 

0.042 

— 

— 

125.00 

— 


60.8 

78.0 


The reduction in thrust from 0.068 Ibf to 0.042 lbf of 38.3% is attributed solely to the tempera- 
ture drop of the system through the following reasoning. Utilizing the well-known formula for 
thrust, F, through a nozzle, the total thrust for our system of two identical nozzles is 


total 


[C f A t P c ] r 


[ C f A . P c 


Jnozzle 2 


( 3 ) 


where 

C f = the thrust coefficient of the nozzles, 

A t = the nozzle throat area, 

and 

P c = the absolute tank chamber pressure. 

Since c f and A t are constant for a given system, thrust varies directly with the absolute tank 
chamber pressure; i.e., 


F 

p~ ~ 2C f A t - constant . 


Therefore, if a thrust, F 1? at a tank temperature, T 1? is known and it is desired to know the thrust, 
f 2 , at any temperature, T 2 , this may be calculated by 
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where the chamber pressures P c2 and P cl are obtained by knowing the temperatures and utilizing 
the equilibrium chart (Figure 2) for the fuel. Therefore, it is seen that as the temperature drops, 
the thrust drops, and vice versa. 

The explanation of the temperature drop and its extent is as follows: Since the reduction in 
balance displacement was approximately linear, by assuming a delivered specific impulse, I , 
of 83 Ibf-sec/lbm, the approximate weight of fuel * f , used, was 


w 


f 




r 


dt 


(0.068 + 0.042) lbf 

r 2(83 lbf -sec/lbm)' ( 18 - 33 sec > = 0.012 1bm. (5) 

The NH 4 HS will attempt to stay on the equilibrium vapor pressure curve shown in Figure 2. When 
0.012 lbm of gas vapor is pulsed out, the solid fuel sublimates to maintain its equilibrium vapor 
pressure. However, the heat of sublimation, h sg , is 782 Btu/lbm and since the fuel tank was 
mounted inside of an 18" x 18" styrofoam container, the process may be considered essentially 
adiabatic. The tank and fuel must provide the sublimation heat; therefore, the fuel temperature, 
and thus equilibrium condition, was necessarily lowered. The heat necessary was 


Q = ( w f)( h J 


= (.012 lbm) (782 Btu/lbm) = 9.38 Btu. (6) 

By assuming that the required heat was obtained from the aluminum tank and NH 4 HS fuel alone, 
the heat balance may be described by: 


Q 



c t dT 


+ 


I’ 

•>0 


w f c f dT 


( 7 ) 


where 

Q = total heat exchanged, 
w t = weight of the aluminum tank (0.610 lb), 
c t = specific heat of aluminum (0.214 Btu/lb. °F), 
w f = weight of NH 4 HS fuel (1.600 lb), 

and 

c £ = specific heat of NH 4 HS (0.340 Btu/lb. °F). 
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RELATIVE DISPLACEMENT RELATIVE DISPLACEMENT 


Assuming constant specific heats and weight of NH 4 HS, and that the fuel and tank equalize at the 
same temperature, the change in temperature should have been approximately 


w t c t+w f c f 

9.380 Btu 

/ Btu \ / Btu \ “ 13.9°F . 

(.610 lb) ( . 214 T&f) + (1.610 lb) (.340 Tb^) 

From Table 2 it is seen that the tank temperature (temp. 1) actually dropped from 73.4°F to 60.8 °F 
or 12.6°F. 




Figure 9— Plot of relative dis- 
placement versus time for a 
250-millisecond pulse bit. 



Figure 10— Plot of relative dis- 
placement versus time for a 
77-mil I isecond pulse bit. 


TIME (seconds) 






Test 2 


The timing circuit external to the chamber supplied 28.0 volts to the solenoid valve for 250- 
millisecond pulses every 3 seconds for a total of 40 pulses. The purpose of this test was to 
simulate the propulsion system incorporated in the pointing controls in a spinning satellite and to 
investigate the reduction in thrust over approximately 2 minutes of operating time. The 3- second 
intervals simulated a satellite spin rate of 20 rpm and the 250-millisecond pulse simulated a con- 
trols system on time of 30° per rotation. The test record of the first six pulses on the EMB 
pendulum is shown in Figure 8. 

To obtain the initial thrust the first swing relative displacement was measured for the first 
pulse, the total impulse delivered read off of the IMB pendulum calibration curve (see Appendix A), 
the impulse adjusted to correct for the extended pulse (see Appendix A), and the impulse divided 
by the pulse bit to give the average thrust. After the 39th pulse, a 20.8-second delay was inter- 
jected. This allowed the IMB pendulum displacement to completely damp out. The 40th pulse was 
then delivered. The relative displacement of the 40th pulse was utilized to indicate the final 
impulse and thrust. Table 3 tabulates the results of Test 2. 

Again the terminal temperatures shown indicate the lowest values of temperature recorded 
and it was assumed equilibrium was obtained at that time. 

From Table 3 it is observed that the average thrust dropped from 0.064 to 0.043, or 32.8%. 
Again, as in Test 1, it is attributed directly to the drop in temperature of the tank from 70.7°F to 
60.8°F. 


Table 3 
Test 2 Data. 


Pulse 

Time 

(seconds) 

Relative 

Displacement 

Total 

Impulse 

(lbf-sec.) 

Thrust, 
2 nozzles 
(ibf) 

Temp 1 
<°F) 

Temp 2 
<°F) 

1 

0 

2.26 

.0159 

.064 

70.7 

82.4 

40 

134.8 

1.52 

.0107 

.043 



— 

300.0 

— 


— 

60.8 

73.5 


Test 3 

The timing circuit external to the vacuum chamber was utilized to supply 28.0 volts to the 
solenoid valve for six impulse bits verying in duration from 250 to 77 milliseconds. The record 
of 250 and 7 7- millisecond pulse bits are shown in Figures 9 and 10. 

The impulses and thrusts were reduced from the IMB calibration curve as described in 
Appendix A. Table 4 summarizes the results of test 3. 
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Table 4 


Test 3 Data. 


Pulse 

Pulse Bit 
Duration 
(milliseconds) 

Relative 

Displacement 

Total 

Impulse 

(lbf-sec.) 

Thrust, 
2 nozzles 
(lbf) 

Temp 1 
(°F) 

Temp 2 
(°F) 

1 

250 

1.50 

0.0106 

0.042 

63.5 

75.2 

2 

200 

1.29 

0.0089 

0.045 

64.4 

75.2 

3 

150 

1.05 

0.0071 

0.048 

65.3 

78.8 

4 

100 

0.80 

0.0054 

0.054 

68.0 

82.4 

5 

77 

0.68 

0.0045 

0.059 

69.8 

84.2 

6 

77 

0.68 

0.0045 

0.059 

69.8 

84.2 


Prior to pulse 1, the heater strip wiring was utilized to elevate the temperature of the fuel 
tank and valve. From Table 4 it is seen that the temperature continued to rise during the six 
pulses; this effected the increase in average thrust as indicated in the table. 

POST TEST EXAMINATION 

After the termination of the tests the propulsion system was left inside the vacuum chamber 
at 2 x 10" 2 mm of Hg for approximately 24 hours. The chamber was then brought back up to 
atmospheric pressure and the propulsion system was removed, purged with dry nitrogen, and 
stored in a small air-tight can at 34°F. 

Two weeks later the fuel system was taken out of storage for examination and the following 
was revealed: 

1. The fuel vapors had attacked the black anodized surface of the aluminum tank causing a 
slight discoloration as may be seen in Figure 1. 

2. The valve exterior surface was nickel plated over mild steel and had been attacked by the 
hydrogen sulfide forming a crust on the metal. 

3. Disassembly of the valve revealed that there were brass components inside the valve 
which had been severely attacked by the hydrogen sulfide such as to render the valve inoperable. 

Since the valve was operating at termination of the tests and none of the exterior effects were 
obvious when the propulsion system was removed from the vacuum chamber, it is surmised that 
the valve developed a leak near the end of the tests. Stored in the small can, the leak allowed high 
concentrations of the fuel vapor to come into contact with the exterior of the system and caused 
the corrosive attack. The existence of a leak, is further supported by the following analysis of the 
total fuel missing versus what should have been used. 
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The total impulse delivered by the fuel system during the three tests was 


Test 1 1.009 lbf-sec 

' Test 2 0.535 lbf-sec 

Test 3 0.041 lbf-sec 

1.585 lbf-sec 


Therefore, the theoretical weight of fuel w f used, was 

2F At 

w f “ “1 

sp 

1.585 lbf-sec 
w f - 83 lbf-sec. /lbm 

w f - 0.019 lbm , 

where 

2FAt = the total impulse delivered, 

and 


I = the specific impulse of the fuel. 

The actual fuel used during the period while the propulsion system was inside the vacuum 
chamber was obtained by weighing the fuel tank before it went in and after it came out. The actual 
total fuel used was 0.197 lbm. The order of magnitude difference in fuel can only be explained by 
a leak. 

CONCLUSIONS 

1. The impulse and thrust derived from the solid sublimating propulsion device did not 
deteriorate in any way due to the extended storage period of 8 months. The thrust data indicated 
that the system met the original design specifications very well. 

2. There is a corrosion problem when high concentrations of ammonium hydrogen sulfide fuel 
are allowed to remain in contact with certain materials for extended periods. This led to the failure 
of the valve. Materials incompatible with this particular fuel, such as brass and copper, could not 
be used with this propulsion system where they may come into contact with the fuel vapors. There 
is some evidence to the effect that the corrosion may be predicated upon the existence of moisture 
or air in which case the problem may be alleviated in outer space. 
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3. The concept of spacecraft control devices utilizing the solid sublimating fuel as a means 
for propulsion is very feasible and attractive. Assuming installation of a proper valve, the atti- 
tude control system utilizing the solid sublimating fuel as the propulsion means is ready for the 
transition to flight hardware. 
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Appendix A 

The Integrating Microthrust Balance 

The Integrating Microthrust Balance, IMB, is a com- 
mercially developed and produced device for accurately 
measuring the jet performance of very low thrust jet reac- 
tion control systems both in the steady state and pulsed 
operating modes. The contractor who developed the balance 
and calibration techniques was Rocket Research Corpora- 
tion. This Appendix will describe the EMB, show its prin- 
ciples and procedures for use, and develop formulas for 
correcting the IMB data output when extended pulses are 
measured with damping in the system. 


13 


General Description 




The IMB basically is a rectangular testing platform supported at the four corners by very 
weak springs such that the platform takes on the characteristics of a pendulum. Figure A1 illus- 
trates the main components. 

The microthrust balance uses a soft suspension to measure the momentum imparted to it by 
a single jet pulse. A linear displacement pendulum concept is employed, which is capable of 

measuring thrust impulse durations from less 
than 10 milliseconds up to 350 milliseconds 
with an overall accuracy of better than 2.0%. 
When the propulsion system to be tested is 
pulsed, the total impulse is absorbed by the 
balance which then deflects in proportion to 
this impulse and then continues to oscillate 
back and forth. The initial amplitude of this 
oscillation is directly proportional to the de- 
livered impulse. The relative displacement is 
therefore measured remotely by an optical 
displacement detector. By calibrating the pendulum balance for relative displacement versus 
known values of impulse, the impulse delivered from a propulsion system may be obtained simply 
by recording its relative displacement and utilizing the momentum versus relative displacement 
calibration curve. 


V////////////////////////^ 



PROPULSION TEST SYSTEM 
SUPPORT ARMS (4) 

SPRING SUSPENSION (8) 

RECTANGULAR 
PLATFORM 


Figure A1 —Side view of integrating microthrust balance. 


Basic Principles 


There are certain assumptions paramount for the ideal operation of this soft balance. They 

are: 


1. The friction and other damping mechanisms acting on the balance are negligible with re- 
spect to the inertial forces, so that, within the desired accuracy, the balance may be treated as 
an undamped system. 

2. The impulse is delivered to the balance rapidly enough so that negligible displacement 
occurs during the thrust period. 

3. The supporting springs have a mass much smaller than the mass of the balance. 

4. The platform and thrust mount located on it have a very high spring constant compared to 
the mounting springs. 

5. The spring constant and mass of the system remain constant during the pulse measurement. 
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If the requirements above are met, the mass- spring system motion is described by the fol- 
lowing simple equations. These are: 


FAt - mv , 


where 

F = motor thrust, 
At = pulse duration, 
m = balance mass, 

and 


v = balance velocity. 

After the pulse is applied the balance has a velocity and therefore a kinetic energy (1/2 mv 2 ). 
This kinetic energy is ultimately all stored as potential energy of the mounting springs at maxi- 
mum deflection according to the relation, 


1/2 mv 2 - 1/2 kx 2 . 


Solving for momentum, 


mv - x(mk) 1/2 , 


where 

k = effective spring constant of the balance pendulum 


and 


x = deflection of the balance pendulum. 

Since m and k remain constant for a given test set-up, it is seen that the balance deflection, x, 
varies directly with the balance momentum, mv. 

Calibration of the Balance 

For the tests described in this paper, the following calibration of the IMB was conducted. 

With the test system on the balance and the entire setup inside a vacuum chamber at 2 x 10' 2 mm 
of mercury, known values of momentum were applied to the balance and the resultant relative 
displacements recorded. The momentum applied was accurately determined by measuring 
the velocity of a ball which was propelled at and captured by the balance pendulum. By 
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Figure A2— Calibration curve showing 
momentum versus displacement. 


dicated impulse may be adjusted or corrected for 
development of the correction curve. 


knowing the ball weight, the applied momen- 
tum, mv, was obtained. Six balls or vacuum 
calibration points were recorded. The cali- 
bration curve for the test is illustrated in 
Figure A2. 

Finite Pulse Bits With a Damped System 

For processes where all five assumptions 
previously mentioned are valid, the impulse 
delivered by any pulse would be that obtained 
from entering the calibration curve at the cor- 
responding displacement. However, the actual 
conditions encountered during this test differed 
from the ideal case in two respects. First, the 
pulse duration and suspension length were such 
that it was impossible to impart all of the 
impulse to the balance before a finite displace- 
ment occurred. Second, the balance pendulum 
motion was damped by allowing a paddle under 
the pendulum to hang in a fixed pot which was 
filled with 20,000 centistoke silicone oil. The 
effect of these two conditions was to slightly 
reduce the balance deflection for any given 
impulse and thus introduce a small error. 
However, by manipulation of the equations of 
motion governing these phenomena, a correc- 
tion curve may be developed by which the in- 
any finite pulse. The following treatment is the 


To develop a correction curve the approach will be to derive the two expressions describing 
the maximum displacement for the two processes involved; that is, first, for the ideal case balance 
displacement, x max , when all of the momentum is imparted to the damped balance prior to any 
displacement and then, second, for the maximum balance displacement, x max , when the balance 
pendulum is displaced and moving while the momentum is being applied. Then for any given pulse 
duration the ratio of x max/ /x majc will indicate the correction necessary to the actual pendulum dis- 
placement to indicate the true total momentum applied. 

The equation of motion of the spring-mass system with damping is 


mx + cx + kx F(t) 


(Al) 
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where 


m = effective mass of the balance pendulum 
c = damping of the balance pendulum 
k = effective spring constant of pendulum 
x = displacement of the balance pendulum 
x = velocity of the balance pendulum 
x = acceleration of the balance pendulum 

and 

F(t) = forcing function dependent on time. 

When all of the momentum is imparted to the balance prior to any displacement, the forcing 
function, F(t), is an impulse function at t equal zero. The solution of the linear differential equa- 
tion for the complimentary function then yields the expression describing the pendulum motion: 

v o 

X - e~ pt — sin qt , (A2) 


where 

t = time, 

V 0 = initial pendulum velocity, x, at t = 0, 

p = c /2m (a substitution made for simplicity throughout this Appendix), 
q = l/2m V 4km - c 1 (a substitution made for simplicity throughout this Appendix). 


By differentiating Equation A2 and setting the result equal to zero, the time, t, for the maximum 
displacement is obtained. Substituting this back into Equation A2, yields an expression for the ideal 
case maximum pendulum displacement, x max , which is 


X 


raaxj 


v o 

Vp 2 + q 2 


-(p/q) sin 

e 


^q/^p 2+ q 2 ) 


(A3) 


V — i 

0 (P/q) sin ( w j/^ n ) 

03 e 

n 


(A4) 


since 


&,• - yp + q 2 = undamped natural frequency of the balance pendulum, 
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and 


oo d = q n damped natural frequency of the balance pendulum. 

Now it is necessary to derive the expression for maximum displacement, x max , when the thrust 
of the rocket, F(t), is spread over a finite time period. 

As the jet pulse width increases, the forcing function F(t) can be approximated with good ac- 
curacy by a pulse of constant amplitude. That is, since the thrust rise and decay transients are 
on the order of a few milliseconds each for pulse widths in the range between 50 and 350 milli- 
seconds, F(t) can be effectively approximated by a single square-wave pulse. 

For pulse widths in the range of 10 to 50 milliseconds, the pulse duration is sufficiently short 
with respect to the balance natural frequency that the assumption of negligible displacement during 
the thrust period is a good one, and the solution given by Equation A4 is valid. 

For extended thrust periods, the general solution to the differential Equation A1 is found to be 


X 


1 - e~ pt |cos qt ~ ~ sin qt^ 


(A5) 


where 


F = thrust of the rocket. 

For an extended thrust period (i.e., F = constant) of time, t ., followed by the thrust going to 
zero (i.e., F - 0), the displacement expression is solved to be 


X 


e~ pt 


i 


v i + pyi 
q 


sin qt ' + Xj cos qt ' 


(A6) 


where 

t : = thrust pulse duration 

t = tj + t ' or t' = the time initiating at instant t l when f - 0, 

Xj = balance pendulum displacement at instant, 1 1 , 

and 

V x = balance pendulum velocity at instant, t x . 

By maximizing Equation A6 the expression for the damped extended pulse displacement is 
obtained; that is, 


X 


max 


“(p/q) s in *D 


V 1 + P X 1 


D + Xj (l-D 2 ) 1/2 


(A7) 
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where 


D = 



1 



(A8) 


It is obvious that expressions for x ; and v t must be developed to evaluate Equation A7. The 
expressions for x s and Vj are derived by utilizing the extended thrust Equation A5 and evaluating 
it at t - tj ; i.e., at t - t,, x ‘ v x> therefore, 


■ ¥ [f ' e "”' 1 ( cos qti + q sinqtj 


(A9) 


and 


v i 



sin qtj . 


(A10) 


We have now derived the two expressions, x max and x max , which describe pendulum displace- 
ment for the two processes in question. Their ratio is 
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(All) 


Inspection of Equations A8, A9, A10, and All reveals that all quantities are either easily measured 
parameters of the balance pendulum or are cancelled out when x 1? v x , and D are substituted into 
Equation A1 1 . 

By inserting all known quantities into the four equations, a correction curve was developed 
which described the finite pulses on our damped system. The curve is shown in Figure A3. In- 
spection of the equations and the correction curve indicates that the error in displacement for 
finite pulse durations chi any given test set-up is directly dependent upon the pulse duration only. 
Therefore, for a given pulse of duration, At a , the correct delivered impulse was obtained by the 
following procedure. The relative balance displacement, S a , was recorded. By entering Figure A2 
at S a , the uncorrected impulse, l unc , on our damped balance was obtained. Then by entering Fig- 
ure A3 at At a , the analytical correction constant, R, for our damped finite pulse was read off. The 
true total impulse, l tot , delivered by the finite pulse was then, 
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tot 


R 


(A12) 
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PULSE DURATION (milliseconds) 


Figure A3— Finite-pulse correction curve showing 
correction constant versus pulse duration. 

Steady State Thrusts 


It should be noted that a much simpler correc- 
tion formula may be obtained if damping is neg- 
lected in the entire analysis. For the damped 
system on which our tests were performed, the 
correction curve developed by neglecting damp- 
ing would be only slightly different from Fig- 
ure A3. If, however, the damping is such that 
there is a significant difference between the 
damped natural frequency, co d9 and undamped 
natural frequency, a> n , then the equations de- 
veloped in this Appendix may be utilized to 
obtain the correct impulse delivered by the 
finite pulse on the damped system. 


An equivalent system to describe steady state thrusts on the pendulum balance will now be 
developed. From the energy balance describing sinusoidal balance motion, that is, 

1/2 mv 2 = 1/2 kx 2 , 

7 max ' max ’ 


the following relationship may readily be developed which correlates the sinusoidal displacement 
resulting from ai* applied momentum to an equivalent displacement which would result from an 
equivalent steady state force system. 


mv 



(A13) 


where 


F = a steady state thrust only 

and 


cc n = undamped natural frequency of the balance. 

Equation A13 indicates that the resultant displacement of the balance from a steady state 
thrust, F, is equal to the maximum balance displacement by an ideally applied momentum (mv 
applied before initial displacement of balance) where the momentum is numerically equal to r/u n . 
Since the balance calibration curve is linear, it is only necessary to know the momentum applied 
and resultant maximum displacement at one point; from this any steady state thrust displacement 
may be utilized to solve for the thrust, F. Thus, the steady state thrust was obtained by the fol- 
lowing procedure. 
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1. Figure A 2 was entered at any displacement and the corresponding momentum read, i.e., 
for Sj = 1.5, mv = 0.01 lbf-sec. 

2. An equivalent steady state thrust, F eq , was calculated which would result in the same dis- 
placement, s x = 1.5, i.e., 

F e , = >nvco n = (0.01 lbf-sec) ( 4.5 
= 0.045 Ibf . 

3. For small displacements, 

F 

— - k - constant , 

and we have 

.045 

k = = 0.03 - 

4. Therefore, with our particular test set-up, for any steady state thrust, F a , and the result- 
ant displacement, S a , the thrust was 

F fl = 0.03 (S.) . 

As an example, if s a equaled 1.40, then 

F a = 0.042 lbf . 
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Appendix B 

List of Symbols 


A t = Throat area of the test rocket 
c = Damping of the balance pendulum, 

C f = Thrust coefficient of the test rocket 
c f = Specific heat of the fuel 
c t = Specific heat of the fuel tank 
F = Thrust from test rocket 
h = Heat of sublimation of the fuel 

s g 

l sp = Specific impulse of the fuel 
k = Thrust coefficient of the test rocket 
m = Thrust balance mass 
P c = Absolute tank chamber pressure 
q = Total heat exchanged 
T = Temperature 

t = Time 
w f = Weight of fuel 
w f = Weight of fuel used per unit time 
w t = Weight of the fuel tank 
x = Thrust balance lateral displacement 
x = Thrust balance lateral velocity 
x = Thrust balance lateral acceleration 
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